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Abstract 



O In the previous paper |arXiv:091 1.06791, we showed that the Reissner-Nordstrom black hole 

O 

^^ in the 5-dimensional anti-de Sitter space coupled to the Maxwell theory with the Chern- 

Simons term is unstable when the Chern-Simons coupling is sufficiently large. In the dual 
conformal field theory, the instability suggests a spatially modulated phase transition. In 
this paper, we construct and analyze non-linear solutions which describe the end-point of 
this phase transition. In the limit where the Chern-Simons coupling is large, we find that 
the phase transition is of the second order with the mean field critical exponent. However, 
the dispersion relation with the Van Hove singularity enhances quantum corrections in the 
bulk, and we argue that this changes the order of the phase transition from the second to 
the first. We compute linear response functions in the non-linear solution and find an infinite 
off-diagonal DC conductivity in the new phase. 



1 Introduction 

In the previous paper [Ij, together with Shin Nakamura, we pointed out that the MaxweU 
theory with the Chern-Simons term in the 5-dimensional Minkowski space is tachyonic when 
a constant electric field is turned on. A similar mechanism leads to an instability of charged 
black hole in the 5-dimensional anti-de Sitter space (AdS^) if the Chern-Simons coupling 
for the Maxwell field is sufficiently large. Interestingly, the instability modes carry non-zero 
momenta along the boundary of AdS^. This suggests that there is a novel phase transition 
in the holographically dual field theory at finite chemical potential, where order parameters 
acquire spatially modulated expectation values. 

The analysis of our previous paper was at the linearized level, and what we observed was 
an onset of the phase transition. To understand the nature of the new phase which emerges as 
the end-point of the instability, we need to examine full non-linear solutions to the equations 
of motion in the bulk. In this paper, we construct such solutions in the limit where the Chern- 
Simons coupling is large and back-reaction of the Maxwell field to the metric is negligible. 
This is analogous to the probe limit employed in [2j. Using the solutions, we compute the 
expectation values of the order parameters near the phase transition temperature and find 
that the phase transition is of the second order with the mean field critical exponent. 

The Chern-Simons term modifies the dispersion relation in such a way that the density of 
states per unit energy diverges at some non-zero momenta, causing the Van Hove singularity. 
Moreover, this happens even above the phase transition temperature. It suggests that quan- 
tum corrections to the phase transition can be significant. We argue that the order of the 
phase transition is changed from the second to the first due to quantum effects in the bulk. 

The rest of the paper is organized as follows. In section 2, as a warm-up exercise, we 
will discuss non-linear solutions to the Maxwell-Chern-Simons theory in the 5-dimensional 
Minkowski space. In section 3, we turn to the theory in the full black hole geometry and 
construct non-linear solutions in the limit where the Chern-Simons coupling is large. We find 
that the phase transition is of the second order with the mean field exponent. In section 



4, we discuss quantum corrections to the phase transition and argue the order of the phase 
transition is changed. We evaluate the hnear response of the system in section 5. In appendix 
A, we discuss non-hnear solutions in AdS2 x K^, which is the near horizon limit of the extremal 
black hole. 

Note added: After the first version of this paper was completed, we were informed of the 
work [3j , which suggested that an instability to crystalline phases might be a generic feature 
of phases which are describable by a bulk AdS2 geometry. Such an instability would provide 
a natural way to understand the ground state entropy. 

2 Maxwell-Chern-Simons Theory in Minkowski Space 

In this section, we consider the Maxwell theory with the Chern-Simons term in the 5- 
dimensional Minkowski space. The Lagrangian is given by 

C = -\fuF'' + ^e'J^^^AiFjKFLM , (2.1) 

where I^ J, . . . run from to 4. The equations of motion are 

djF'i + ^e^'^'^^'Fj^FLM = . (2.2) 

In particular, the time component of the above can be written as the Gauss law, 

a^n^ = , where H^ = -F^° + ae^^^^'AsFcD , (2.3) 

where the indices A, S, . . . run from 1 to 4. 

A constant electric field is a solution to the equations of motion. However, as shown in 
[T], there are unstable modes in the following range of momentum, 

< |fc| < 4aE, (2.4) 

where E is the background electric field and fc is a projection of the spatial momentum onto 
the plane orthogonal the electric field. Let us describe the instability mode found in [Ij. If the 
electric field is in the x^ direction, it is convenient to decompose the 5-dimensional momentum 



as (p^, fc^), /x = 0, 1 and i = 2, 3, 4. Consider a linear fluctuation of the Maxwell fleld of the 
form, 

a, = cJ^^e^^^"^+^^^"\ (2.5) 

where c\ ^ are an eigenvectors of eijkkj^ 

e,,,k,cf^ = ±\k\cl^\ (2.6) 



obeying the transverse gauge condition kiC- = 0. Substituting (2.5) into the equations of 



motion (2.2), we flnd the dispersion relation for this mode as 

bo)' - {Pif = k'T^aEk , k=\k\. (2.7) 



This means that it is tachyonic for the range (2.4) if we choose c- . 

2.1 Non-linear Solutions 

We can flnd a non-linear solution triggered by the perturbation a^ = q ^e^^^^^+^^i^' as follows. 
Although the unstable mode breaks the translational invariance along the direction of the 
momentum fc, it is invariant under some combination of translation along k and rotation in 
the transverse plane. The solution is also translationally invariant in the transverse directions. 
It is then natural to look for a non-linear solution with the same set of symmetries, and we 
choose the following ansatz, 

^ = , A2 = g{t) 

Ai = A,{t) (2.8) 

As + iA^ = h{t)e-'^''' . 
We denote the time coordinate by x^ or t interchangeably. The equation of motion for the 
function g{t) sets g{t) = 0. The remaining equations of motion become 

Ai{t) -Aakh{t)h{t) =0 

(2.9) 
h{t) + k^h{t) + AakAi{t)h{t) = . 

Note that the momentum Hi conjugate to Ai is given by 

Hi = -ii(t) + 2akh{tf . (2.10) 



The first equation of (2.9) can be written as dUi/dt = and solved by Hi = constant. 



The integration constant is fixed as Hi = E by the initial configuration where h = and 
$1 = -Ai = E. 



From (2.9) and (2.10), 



h{t) + r/i(t) - 4akEh{t) + 8a^k^h\t) = . 



(2.11) 



This can be integrated and yields 

-h{tf - -k{4:aE - k)h{tf + 2a^k'^h{tf = constant. 

Zi Zu 

This is equal to the energy density ?{ of the electro-magnetic field, 

'H = -F^aF^a + jFabFab 5 



(2.12) 



(2.13) 



minus the energy density ^E'^ of the constant electric field. Thus, we are eflFectively considering 



a classical particle with coordinate h{t) moving in the potential 

U = --k{AaE - k)h^ + 2a^k^h^ . 

Zu 



(2.14) 



This is a double well potential for < A; < ^aE as in Figure [T] The original homogeneous 
phase corresponds to the point /i = 0, which is unstable. If we add some perturbation, the 
amplitude h{i) starts oscillating as in the figure. 




Figure 1: Double well potential U for a classical particle with coordinate h. If the particle 
starts slightly outside of the origin, say at ^4, then the particle will oscillate between A and 
another point B with the same potential energy. 



2.2 The Final Configuration 



We have seen that the instabihty induces an osciUatory solution in the potential (2.14). 



Suppose that our system is weakly coupled to a heat reservoir with a large heat capacity at 

very low temperature. Eventually the oscillation will fade away by transferring its energy to 

the heat reservoir and the system will land on its lowest energy state. Let us try to find out 

the final configuration of this process. 

For the static solutions, 

, 1 UaE - k ,^_, 

which stay at one of the two minima of the potential, the energy density H is given by 



1 ^ k f k 

2 4a V 8a 



(2.16) 



Especially, the energy density vanishes when k = 0. Figure [2] shows the amplitude h and 
the energy density H change as functions of k. Since the energy density is monotonically 
increasing in fc, we expect that solutions with k y^ are unstable and decay to the solution 
with k = 0. Note that, although h diverges as k goes to 0, the field strength E vanishes in 
this limit. The constant electric field in the initial configuration is wiped out in the k = 
solution and the final configuration will be the trivial vacuum state with E = 0. 





4aE 



4aE 



Figure 2: The amplitude h and the energy density ?{ as a function of k. 



We can directly check that the static solution (2.15) with fc 7^ is unstable. The solution 



i(o)/ 



l(0). 



l(0) 



corresponds to the gauge field configuration, A\ \t) = ^t and A\ +iA\ ^ = ^\/^^^-^e ^^^^ 



Let us add a small perturbation A^^^ -^ A^^^ +a to this background. Assuming that the modes 
depend only on t and X2 with ao = 0, the equations of motion become 

Dai + 2akh{e'^'''doa+ + e-'^'^'doa^) = 

Das = (2.17) 

na+ + ikd2a^ - 4:akhe~'^''^doai = , 

where D = — Sq +5| and a+ = a^ + ia^. The coefficients of the equations have X2 dependence, 
which can be removed by using two real variables a^ and 0.4 such that a+ = e~^^^'^{as + 10.4) 
in place of a. Assuming the t and X2 dependence of the fields ai, ds and a^ to be of the form 
^-iujt+iqx2 ^^^ |-]^g^^ (22 = 0, a non-trivial solution exists if and only if 

(cj2 - q^) [cj^ - 2uj\q^ + fc(4a£; - A;)) + q\q^ - fc^)] = , (2.18) 

From the second factor, we see that the product of two solutions for cj^ is q^{q^ — fc^), which 
is negative for < q < k. Thus one of the two roots of u'^ must be negative, representing an 
unstable mode. Since the momentum q of the instability mode is in the range < g < fc, we 
expect that the solution decays toward the lowest energy state with momentum k = 0. 

We have also performed numerical analysis of time dependent solutions with the initial 
configuration of constant electric field E assuming that solutions depend only on the coor- 
dinates t and X2- We found that a small localized perturbation generates a domain where 
the electric and magnetic fields fiuctuate and that the domain expands at the speed of light. 
The magnetic field in the domain carries a range of momenta, which tend to move to zero 
momentum state. The strength of the electric field decays as the domain expands, suggesting 
that the system will eventually settle down to the trivial state with E = 0. 

To summarize, the instability of the constant electric field in the Maxwell-Chern-Simons 
theory in the 5-dimensional Minkowski space, which we found in our previous paper [Ij, leads 
to the trivial vacuum state with no background field strength Fjj = 0. This reminds us of 
the Schwinger mechanism where a constant electric field is screened by virtual production of 
electron-positron pairs. 



This result should be contrasted with the corresponding instability of the charged black 
hole in AdS^^ which we will study in the next section. There, we do not expect that the 
background electric field to disappear since the electric charge of the black hole is fixed by 
the chemical potential at the boundary. Indeed, we will find stable solutions with non-zero 
momentum k in this case. 

3 Maxwell-Chern-Simons Theory in the AdS^ Black Hole 
Geometry 

In this section, we will construct non-linear solutions which describe the end-point of the 
instability of the charged black hole in AdS^. Since the Maxwell field contributes to the 
energy momentum tensor, in general we need to analyze the coupled Einstein and Maxwell 
equations. Here we will simplify the problem by taking a limit where we can ignore the 
back-reaction of the Maxwell field to the metric. 

The Lagrangian density for the Maxwell-Chern-Simons theory is given by 

C = -\fjjF'' + .^^^i'^^^MFjKF^M . (3.1) 

Rescaling the gauge field as yl = aA^ the Lagrangian density becomes 

/: = A f-7^^^^" + i^^^'^'^^MFjkF^^ . (3.2) 



When a is large, for a solution with finite A^ the energy-momentum tensor is of the order 
0{\jo?^ and the coupling of the Maxwell field to the metric is suppressed. This limit is analo- 
gous to the infinite charge limit considered in the holographic description of superconductivity 

12J. 

3.1 The Large ol Limit 

To keep A finite, we have to scale the background gauge field as well. This means that the 
chemical potential /x of the black hole should also be scaled in such a way that the combination 
a\i remains finite. Let us examine what this limit means to the black hole solution. The 
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Reissner-Nordstrom solution has the metric 



ds^ = -H(r)dt^ + ^r---dr^ + r^dx^ , x = (x^ x^ x^) , (3.3) 

Hir) 



where the function H{r) is given by 



H{r) = r^ 
The temperature in this hmit of /i ^ becomes 



-r-±) 

r J 



(3.4) 



^+ /o MM . ^+ 



r = -^ 2 - f^ ^ ^ . (3.5) 

The background geometry in this hmit is simply the (neutral) Schwarzschild AdS^ solution. 
In terms of the rescaled finite gauge field A^ the background field strength is given by 



—dt A dr , where E = —2aiJir^ = ^- 

J.6 -T ^ ^ 



^dt Adr ^ where E = —2aiJir]^ = — Z~Z^ (^•^) 



where we introduce a new variable r = rj^/Tiaji for later convenience. Since r = Tj jia^ it can 
be thought of as a rescaled temperature. 

It is important to note that we have access to the phase transition point in this limit. In 
[1], we studied the instability of the Reissner-Nordstrom solution and obtained the critical 
temperature T^ for the instability as a function of the Chern-Simons coupling a. The result of 
our numerical analysis is reproduced in Figure |3l For large a, the dimensionless combination 
Tc/ jJi grows linearly in a. Thus, we can analyze the behavior of the system near T = Tc by 
taking the limit of a ^ oc while keeping the combination aji finite. 

Let us find a non-linear solution to the equations of motion that describes the spatially 
modulated phase in this limit. We look for a solution that has the same symmetry as that of 
the unstable modes found in [Ij, namely a linear combination of a translation along X2 and a 
rotation in the 3-4 plane, as well as the translation symmetries along t, X3 and X4. This leads 
to the following ansatz, 

^ = /(r) , ii=^i(r) , M = g2{r) 

. . . , (3.7) 

A'.= A^ + iA^ = h{r)e-'^''' . 
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T ^ Tc{a) 




0.05 - 



fl/^crit 



Figure 3: Critical temperature as a function of the Chern-Simons coupling a. The shaded 
region indicates a phase with a non-zero expectation value of the conserved current J which 
is helical and position dependent. 

Note that gi{r) can be set to vanish by a gauge choice, and 5^2 (^) has to vanish by the Maxwell 
equation jj- = 0. The non-trivial equations of motion are 



-drir'^fir)) - 4kh{r)h\r) = 

dr(H(r)rh\r)) - -k^h(r) + ikh(r)f(r) = . 
r 



The first equation can be integrated and becomes 



(3.8) 



rV (r) + 2kh{ry = E 



(3.9) 



Eliminating f'{r) in the second equation by using the above relation, we obtain 



1 



h{r) 



dr{H{r)rh'{r)) - -k^h{r) + 4k^{E - 2kh{rY) = 



3.2 Second Order Phase Transition 



(3.10) 



We have solved the differential equation (3.10) numerically. For each initial condition at the 



horizon, the equation is integrated numerically toward the AdS^ boundary. In general, we 
find a linear combination of normalizable and non-normalizable modes near the boundary. 
Since the new phase of the system should be represented by a normalizable solution, we tune 
the initial condition at the horizon so that the non-normalizable component vanishes. Figure 
^describes numerical solutions for r = 0.35. The left graph shows the amplitude h{r^) = h^ 
at the horizon as a function of the momentum k. 



hlr^ 



0.15 

0.1 

0.05 




A^ 



-0.002 



-*- -k/r+ _o.004 




-k/r+ 



12 3 4 

Figure 4: The amplitude /iq and the energy density difference A'H from the homogeneous 
phase as functions of k. 

Since there is a family of solutions parametrized by the momentum /c, we need to choose 
the minimum energy density solution as the final state. The energy density is given as a sum 
of the electric and magnetic energy 



n 



1 



2^2 



dr 



f\rf + ^-h{rf + 



H{r) ,, 



h\ry 



(3.11) 



It is convenient to rescale the energy density as ?{ = a^r\l-L^ which is finite in the limit of 
a ^ oo. The right graph of Figure |4] shows the energy density as a function of k. Choosing 
the momentum corresponding to the minimum of the energy density, the expectation value 
of the order parameter \{J^ + iJA)\ can be read off from the asymptotic behavior of the 
corresponding bulk field h{r). Figure [5] shows the expectation value | (J3 + iJ/^) \ as a function 
of the rescaled temperature r. Near the critical temperature, it behaves as 



(J3 + ^^4)|=^ 1 



T 



(3.12) 



where A = 2.04 and Tc = 0.37. The critical exponent 1/2 is typical for a mean field theory. 
Indeed, this can be expected from the absence of quadratic terms in the equations of motion 



(3.10) and the fact that we consider the gravity system classically. The mean field behavior 



is also observed in the holographic description of superconductivity [2^ 4J. 

4 Quantum Corrections 

We found that the phase transition is of the second order in the classical supergravity approx- 
imation. In this section, we provide evidence that quantum corrections in the bulk change it 
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|<J3 +iJ4>| 




0.32 



0.34 



0.36 



Figure 5: The expectation value of the order parameter as a function of the temperature. 
The dotted curve is the numerical result and the solid curve is its fit with (1 — r jr^^l'^ . 

to the first order. Such a phenomenon has been observed by Brazovskii [5j and elaborated by 
Swift and Hohenberg in [6j, in the context of a classical statistical model at finite temperature. 
We will extend this result to the gravity theory in AdS^. 

Let us review the Brazovskii model. It is a classical field theory in d space dimensions at 
finite temperature with the following scalar field Hamiltonian in the momentum representa- 
tion, 



T ^ 



d^q 
2.1 {27rY 

A f d^qi d^q2 d^qs 



+ 



0(gi)0(g2)0(g3)0(-gi - 92 - qs) . 



(4.1) 



4! J (2^)^ (2^)^ {27Ty 

Note the unconventional dispersion relation u) = m'^ + {\q\ — qo^. It has the Van Hove 
singularity duj/dq = Sit \q\ = qo^ which will play an important role in the following. 

In the mean field approximation, the system undergoes a phase transition at m^ = 0, and 
there is a spatially modulated phase for m^ < 0. The phase is characterized by the position 
dependent expectation value of (/>(x), 



((/)(x)) = 2a cos (f- X 



(4.2) 



with \q\ = qo- The phase transition is of the second order in this approximation. 

Let us examine if this picture is modified by thermal fiuctuations. The inverse suscepti- 
bility is defined by 

M2 := G-\\q] = qo) = m' - E(|g1 = qo) , (4.3) 



11 



where E represents thermal loop contributions. In the Hartree approximation, 



(4.4) 



where s = irSaq^ '^/(27r)^ with Sa the area of the (d — l)-sphere. There are higher order 
corrections, but they wiU not affect the behavior of M^ near a = 0. The first term in the 



right-hand side of (4.4) is the tree level value. The second term comes from a loop diagram 



as in Figure [6} which gives a contribution near q = q^ oi the form 

A ( d<^Q A ,_, r- dq ^ A TvS.qy ,,.. 

{2nr J M' + {\q\ - qoy ^ {2nr '^' Jo M^ + {q - q,f {2i:Y M ' ^"^ 

There are subleading terms for small M, but what is relevant for the analysis below is the 
1/M pole from the loop. The third term Xa^ comes from the quartic coupling combined with 



the expectation value (4.2) of (f>. 




Figure 6: A loop diagram that contributes to the two point correlation function. 



The free energy J^ for the expectation value (4.2) can be evaluated by setting [S], 



dT 
da 



= 2M^a - Xa^ . 



Combining this with 



da M 



dM Xa 



sX 
2M3 



derived from (4.4), we find 



— f dJ^ , f dJ^ da „ ^ 

J-= / -—da= / - — -—dM 
J da J da dM 

= 2aU + "^^+3^^^-^-2M^) 



(4.6) 



(4.7) 



(4.8) 
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Figure 7: Different shapes of the graphs of the free energy J^ as a function of a as the parameter 
vp? changes. 



Note that the terms in the parenthesis in (4.8) depend on s and A only in the combination 
s\. If we rescale vn and M by a factor of e and s\ by e^, the terms in the parenthesis in 
(4.8) scales like e^. Thus, we can set 5 = A = 1 without loss of generality if we do not care 



about an overall factor of T . Figure [7| illustrates the free energy J^ as a function of m^. For 
m^ > 0, we have a graph with only a minimum at the origin. As w? decreases and crosses 

2 

— 2~33 ^ —1.89, we see two inflection points away from the origin, representing the spatially 
modulated phase. At this point, the homogeneous phase at a = still has a lower energy. As 
we lower vr? further, the energy of the spatially modulated phase becomes lower than that of 
the homogeneous phase, and the flrst order phase transition occurs. 

The crucial point is that the origin a = remains semi-stable throughout the process since 



(Pj^/da^ = 2M^ > at a = 0. This in turn is due to the fact that (4.4) has a solution with 
M^ > for any value of m^ including m^ < 0. This is possible since the second term in the 
equation is proportional to 1/M and diverges as M ^ 0. This is illustrated in Figure [8| Since 
a = is always semi-stable, the phase transition cannot be of the second order. 

The 1/M singularity in ( |4.4[ ) originates from the fact that the kinetic operator [M^ + {\q\ — QoY 
vanishes at M^ = and \q\ = Qq. If the zero of the kinetic operator were at g = 0, it would 
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Figure 8: The relation between m^ and M^. Note that for any value of m^, there is always a 
solution with M^ > 0. 

not have caused the 1/M singularity because the factor q^~^ in the phase space volume q^~^dq 
would have suppressed it. The singularity is generated in our case because of the larger phase 



space volume ^ Qq ^dq at \q\ = qo in (4.5). It ensures that there is a positive M^ solution to 



(4.4) and that the origin of the field space is meta-stable. 



This feature is shared by the gravity theory considered in this paper. In our classical 
analysis in the last section, we found that the phase transition happens at T = Tc and the 
new phase for T < Tc is represented by the non-linear static solutions constructed in section 
[3j Since the phase transition is of the second order, the size of the non-linear solution grows 



linearly in \[T^^^ . This means that, dit T = Tc^ there are static solutions to the linearized 
equations of motion, namely a zero of the kinetic operator at non-zero momenta. As in the 
case of the Brazovskii model, it generates a 1/M singularity in the two-point function. Thus, 
we expect that the homogeneous phase will remain meta-stable for T < T^. If quantum 
corrections are parametrically suppressed (e.g., by 1/A^), the spatially modulated phase will 
eventually acquire lower energy and the first order phase transition will take place at that 
point. 

5 Linear response 

Finally, let us examine linear response of our solution when we couple a gauge field to the 
current J at the boundary. Because of the Chern-Simons term, the current is anomalous. We 
therefore treat the boundary gauge field as an external and non-dynamical source as in recent 
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papers, for example [TIE]. 

Note that the background solution we consider is inhomogeneous, carrying the momentum 
k. In the homogeneous setup, it is natural to choose a translationally invariant source. In our 
case, we consider a small perturbation of the form A = A^^"^ + a, where A^^"^ is our non-linear 
solution and a is a small perturbation with non-zero components, 

as{x2,r,t) + ia4{x2,r,t) = -iaT{r,t)e~'^^^ , a2(r, t) = aL{r,t) . (5.1) 

Here, aT{r^t) and aL{r^t) are real functions of r and t. Notice that we turn on a magnetic 
field as well as an electric field. It is not possible to turn on only an electric field due to 
the Bianchi identity 9[+F^2] = 0, where + denotes the coordinate X3 + ix4. However, in the 
current setup, the magnetic field is determined by the electric field and is not an independent 
quantity. 

The non-trivial linear equations of motion of the fields aT(r, t) and aL{r^ t) are 



r 



H{r) 



diaL + dr[rH{r)draL) + idrh{r)dtaT = 



^ o.,uuf.^2 (5.2) 



"^ -d^ar + dr{rH{r)draT) - —ar + 4fc ^ ^^kh{r) ^^ _ ^q^^^^^q^^^ ^ q 



H{r) r r^ 

Here h{r) is the non-linear solution we have found previously. We are interested in modes 
with definite frequency of the form e~^^^ The fields aT{r) and aL{r) behave as H{r)~^ near 
the horizon r = r+ if we impose the in-going boundary condition. On the other hand, near 
the AdS^ boundary r = oc, they behave as 

ar = aT(o) + aT(2)r~^ + <^T(iog)^"^ log r/r+ H 

(5.3) 
ctL = ciL(o) + ciL(2)r ^ + aL(\og)r ^ log r/r+ H . 

Note that there are logarithmic terms and that the coefficient aT(2) or aL{2) can be shifted by 
<^T(iog) and aL(iog) if we change the scale of the radial coordinate r. This corresponds to the 
choice of a renormalization scale, and should not affect physical quantities. The behavior of 



the equations (5.2) near the AdS^ boundary r = oc fixes ^^(log) and ^^(log) to be 



1 CJ^ 

<^T(log) = 2*^^ ~ ^ )<^T(0) , ^L(log) = ^<^L(0) • (5.4) 
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Since we are solving a set of linear differential equations, the coefficients aT(2) and ai,(2) are 
determined linearly from aT(o) and a^oy That is, 

\aL(2)J \amj 

for some 2x2 complex matrix fl. From the prescription in [U], and analogously to [TU], the 
retarded Green's functions for fields gt and gl are given by 




,/ °. .N h (5-6) 



.2 |^r.r.JL - I 



where the last term is due to the logarithmically divergent terms in (5.3) whose coefficients 



are given in (5.4). These divergent parts can be removed by adding a suitable boundary 



counterterm to the gravity action. After subtracting the divergence, we arrive at 

G^ = 2n- (^2 (^' - ^') ,°,) , (5.7) 

When considering modes of definite frequency, the additional electric fields Et and E^ are 
i(jjaT(o) and iuja^i^Q)^ respectively. Therefore the conductivity is given by 

Let us evaluate the conductivity a numerically for the Chern-Simons coupling a = 0.59. 
For this value of a, the minimum energy occurs at the momentum k = — 2.37r+. We denote 
the components of the retarded Green's function and the conductivity respectively as 

yR _ I Gtt Gtl 



and 



^^Hs:; s::)' (^-9) 



a = h^ '"^^l . (5.10) 

This numerical result has several interesting features. The imaginary part of the retarded 
Green's function Gtt for a pair of a^ fields behaves like l/u near cu = 0, as shown in Figure 
|9) The pole at w = is directly related to the fact that a^ is a Goldstone mode. This can 
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Figure 9: The real and the imaginary parts of G^^ 
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Figure 10: The real and the imaginary parts of aTL and aLT- 



be seen as follows. As u) goes to 0, the solution becomes static and we can ignore the terms 



with time derivatives in (5.2). In this limit, the equation of motion for qt reduces to the 



linearized version of the equation for the phase rotation of the inhomogeneous background 



solution (3.7). Especially, a^ goes to h{r) as we take a static limit. Since h{r) is a solution 
that does not have a non-normalizable mode, if we start with a source such that aT(o) = 0, 



(5.5) implies that Qtt should diverge in the static limit. Therefore, we expect to see a pole 



at cj = in G^^. 



Figure 10 shows the off-diagonal components of the conductivity Im cr^L and Ima^T- Note 
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that they behave as \juj with opposite coefficients. The antisymmetry of gtl and ct^t is due 
to the fact that ( |5.2| ) remains invariant under the change t ^ — t if we simultaneously change 
the sign of aT{r). The poles at cj = in the imaginary part of the conductivity indicate that 
there are delta function contributions at cj = to Rea^L and ReaLT^ which implies that 
there is an off-diagonal infinite DC conductivity. 

Re (diagonal cr) 




. m/r+ 



Figure 11: The real parts of the diagonal components of the diagonalized conductivity matrix. 



It is also instructive to diagonalize the conductivity matrix. Figure [IT] shows the real 
parts of the two diagonal components of the diagonalized conductivity matrix. The real part 
of a conductivity is directly related to the spectral density, and the real parts of the two 
eigenvalues measure the spectral densities corresponding to two linear combinations of the 
transverse and longitudinal currents associated with a^ and a^ fields. In particular it should 



always be positive, and indeed this is shown explicitly in Figure [TT] One of the diagonal 
conductivities has a l/cj^ pole, which is again related to the presence of a Goldstone mode. 
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A Analysis in the AdS2 x ^^ geometry 

In this appendix, we consider the Maxwell field in the AdS2 x K^ geometry ignoring its back- 
reaction to the metric. Since the geometry corresponds to the zero temperature limit of the 
black hole, it does not appear in the limit we consider in section 3. On the other hand, the 
equations of motion can be solved analytically in this case, and it gives a natural generalization 
of the analysis in section 2. 

The metric of AdS2 x K^ is given by 

dv 
ds^ = -r^dt^ + ^ + dxl + dxj + dxl . (A.l) 

Note that, especially, ^/^ = 1. We assume the radius of curvature for AdS2 is 1. A different 
value of the radius can be easily considered by using dimensional analysis. For a generic 
metric, the Lagrangian is given by 

C = -]f:jF'-^ + -^e^'^^'^AjFjKF^^M , (A.2) 

whose equation of motion is 

dAV^gF'') + le^'^^^^F.^Fz.M = . (A.3) 

The index is such that t = and r = 1. Here ^/^ = 1. 

Let us find out a class of non-linear static solutions with the ansatz 

Ao = fir) , Ai = A2 = 
^ := ^3 + iA4 = /i(r)e-^'=^2 _ 
The equations of motion can be written as 

-/"(r) - Aakh{r)h'{r) = 

dr{r^h'{r)) - k'^h{r) + 4akh{r)f'{r) = . 
19 



(A.4) 



(A.5) 



Integrating the first equation, we obtain 

/(r) + 2akh{rf = E (A.6) 

where E is the constant electric background field since the left hand side is just the first 



component of the conjugate momentum Hi. Solving for f\r) and plugging into (A. 5), 

dr{r^h\r)) - k^h{r) + Aakh[r)[E - 2akh{rf) = . (A.7) 

For u = logr, the equation becomes 

h%u) + h\u) - k^h{u) + 4.akh{u){E - 2akh{uf) = , (A.8) 

If we treat 2i as a time coordinate, this equation describes a one-dimensional motion of a 
particle parametrized by h subject to a frictional force and under the potential 

U = -k{iaE - k)h^ - 2a^k^h^ . (A.9) 

If < A; < AaE^ the potential is an upside-down Mexican hat. If a particle starts at one of 
the two hills at 2/ = — oc, it will oscillate around /i = 0, as in the case of the Minkowski space 
discussed in section 2. The new feature in the AdS2 case is that there is a friction term in 
(A.8) and the motion will eventually stop at /i = 0. 

In order for a non-linear solution to exist, the momentum k must obey the additional 
condition, k{k — AaE) > — |, which is equivalent to the Breitenlohner-Freedman bound. This 
condition is needed since any non-linear solution h tends to /i = for large u = log r and obeys 
the linearized equation near the boundary of AdS2- This in particular means that there is no 
non-linear solution with k = 0. Unlike the case of the Minkowski space, the end-point of the 
instability is not a trivial vacuum but a non-trivial solution carrying a non-zero momentum. 
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